Utv ^ 09/787212 

^ ^ / r * , A JCO^'d PCT/PTO 1 4 MAR 2001 

DESCRIPTION 
METHOD FOR FABRICATING A FRACTAL STRUCTURE 

Technical Field 

This invention relates to a method for fabricating 
a fractal structure, especially suitable for application 
to constructing a complex network such as neural network, 
for example. 

Background Art 

Conventional Neuman computers have exercised their 
great power in computers configured to sequentially execute 
specific algorithms and have supported modern scientific 
technologies. In order to design CPU of such a computer, 
a complicated electronic circuit has to be produced and 
progressively optimized by executing its simulation. As 
its design tool, a CAD system for making a complicated 
electronic circuit is indispensable. 

In recent years, information processing learning 
from brains , such as neural network model , has been widely 
researched with a hope of realization as exercising its power 
in pattern recognition. In case a neural network model is 
practically applied in form of a device, it is preferable 
to realize a network such as neural network of a brain in 
form of a certain system . Following such a plan , experiments 
are being conducted toward artificially making nerves for 
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living bodies, and their future development is being 
expected. 

Nerve cells of brains individually have 
complicated tree-like projections to form a fractal 
structure. Such fractal elements grow while interacting 
with each other, and make up a complicated brain neural 
network. 

In order to simulate the function of a preferable 
network , a technique for creating such a brain neural network 
is indispensable. That is, there is a strong demand for 
a technique for making a complicated network, which 
corresponds to a CAD system that has been indispensable for 
conventional Neuman computers. However, there are no 
conventional techniques that make up a structure coupling 
a plurality of fractal elements while controlling their 
fractal nature. 

It is therefore an object of the invention to provide 
a method for making a fractal structure, which can make a 
complicated network like a neural network easily in a 
well-controlled manner. 

A further object of the invention is to provide 
a method for making a fractal structure, which can control 
the coupling mode among different fractal structures and 
can make complicated networks with a more variety of 
structures such as neural networks easily in a 
well-controlled manner. 



Disclosure of Invention 

The Inventors made active studies toward solution 
of the above-mentioned problem, and as a result , found that 
it would be possible to make complicated networks such as 
neural networks by growing fractal structures from different 
start sites and have them interact with each other. 
Additionally, the Inventors found that it would be possible 
to control the coupling mode among fractal structures by 
introducing anisotropy into spaces for growing the fractal 
structures. Thus the Inventors have reached the present 
invention. 

In order to solve the above mentioned subject, 
according to the invention, there is provided a method for 
fabricating a fractal structure characterized in growing 
fractal structures from a plurality of start sites, 
respectively, while having said fractal structures interact 
with each other, to form fractal structures coupled to each 
other. 

In the present invention, universal interaction 
among element fractal structures is controlled against 
fluctuation in growth process of fractal structures from 
individual start sites. 

In the present invention, growth rate from a 
specific start site among a plurality of start sites is 
determined by the probability that a material reaches a 
portion already grown from a remote position due to a 
diffusion process and the probability that a grov/th 



promotion factor reaches the portion already grown from 
portions grown from the other start sites due to a diffusion 
process. The growth rate is proportional to the product 
of the power function of the probability that the material 
reaches the portion already grown from a remote position 
due to a diffusion process and the power function of the 
probability that the growth promotion factor reaches the 
portion already grown from portions grown from the other 
start sites due to a diffusion process* Further, in the 
present invention, fractal property of a structure, its 
self -similarity , complexity or number of coupling is 
typically controllable parametrically. More 
specifically, relative potential determining diffusion of 
a growth promotion factor among individual fractal 
structures grown from a plurality of start sites is 
appropriately adjusted relative to a site at infinity, and 
thereby, fractal property, self -similarity , complexity or 
number of coupling of a structure is controllable 
substantially parametrically. 

In the present invention, anisotropy may be . 
introduced into the space in which the fractal structure 
is grown . More specifically , for example , in case the growth 
rate from a specific start site among a plurality of start 
sites is determined by the probability that a material 
reaches a portion already grown from a remote position due 
to a diffusion process and the probability that a growth 
promotion factor reaches the portion already grown from 



portions grown from the other start sites due to a diffusion 
process, diffusion coefficient in the diffusion process in 
the space for growing the fractal structure has an 
anisotropy. Also in this case, fractal property of a 
structure, its self -similarity , complexity or number of 
coupling is typically controllable parametrically . The 
anisotropy is not limited to the anisotropy of the diffusion 
coefficient, but may be, for example, an anisotropy of the 
dielectric constant in the space for growing the fractal 
structure . 

According to the invention having the 
configuration summarized above, by growing fractal 
structures from a plurality of start sites while having 
individual structures frown from their respective start 
sites to interact with each other to form fractal structures 
coupled to each other, it is possible to fabricate a network 
made up of elements individually having fractal complexity 
easily in a well-controlled manner. 

Further, in case an anisotropy is introduced into 
the space for growing the fractal structure in, when 
individual fractal structures grow from their respective 
start sites while interacting with each other to form fractal 
structures coupled to each other, the response produced 
inside a single fractal structure is affected by a stronger 
nonlinear effect from coupling of different fractal 
structures. Then, remarking the number of sites where 
diff erent fractal structures contact as an index of control 



of the coupling mode among fractal structures, nonlineairity 
can be controlled by controlling the number of those sites . 

Brief Description of Drawings 

Fig, 1 is a schematic diagram that shows a result 
of simulation for fabricating a neural network according 
to the invention; Fig. 2 is a schematic diagram that shows 
a result of simulation for fabricating a neural network 
according to the invention; Fig. 3 is a schematic diagram 
that shows a result of simulation for fabricating a neural 
network according to the invention; Fig. 4 is a schematic 
diagram that shows a result of simulation for fabricating 
a neural network according to the invention; Fig. 5 is a 
schematic diagram that shows a result of simulation for 
fabricating a neural network according to the invention; 
Fig. 6 is a schematic diagram that shows a result of 
simulation for fabricating a neural network according to 
the invention; Fig. 7 is a schematic diagram that shows a 
result of simulation for fabricating a neural network 
according to the invention; Fig. 8 is a schematic diagram 
that shows a result of simulation for fabricating a neural 
network according to the invention; Fijg. 9 is a schematic 
diagram that shows a result of simulation for fabricating 
a neural network according to the invention; Fig. 10 is a 
schematic diagram that shows a result of simulation for 
fabricating a neural network according to the invention; 
Fig. 11 is a schematic diagram that shows a result of 



simulation for fabricating a neural network according to 
the invention; Fig, 12 is a schematic diagram that shows 
a result of simulation for fabricating a neural network 
according to the invention; Fig. 13 is a schematic diagram 
5 that shows a result of simulation for fabricating a neural 

network according to the invention; Fig. 14 is a schematic 
diagram that shows a result of simulation for fabricating 
a neural network according to the invention; Fig. 15 is a 
schematic diagram that shows a result of simulation for 
10 fabricating a neural network according to the invention; 

Fig. 16 is a schematic diagram that shows a result of 
if| simulation for fabricating a neural network according to 

M the invention; Fig. 17 is a schematic diagram that shows 

«| a result of simulation for fabricating a neural network 

ajl according to the invention; Fig. 18 is a schematic diagram 

L|- that shows a result of simulation for fabricating a neural 

"I network according to the invention; Fig, 19 is a schematic 

M 

Cft diagram that shows a result of simulation for fabricating 

q . ■ 

M a neural network according to the invention; Fig. 20 xs a 

20 schematic diagram that shows a result of simulation for 

fabricating a neural network according to the invention; 
Fig. 21 is a schematic diagram that shows a result of 
simulation for fabricating a neural network according to 
the invention; Fig. 22 is a schematic diagram that shows 
25 a result of simulation for fabricating a neural network 

according to the invention; Fig. 23 is a schematic diagram 
that shows a result of simulation for fabricating a neural 
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network according to the invention; Fig, 24 is a schematic 
diagram that shows a result of simulation for fabricating 
a neural network according to the invention; Fig. 25 is a 
schematic diagram that shows a result of simulation for 
fabricating a neural network according to the invention; 
Fig. 26 is a schematic diagram that shows a result of 
simulation for fabricating a neural network according to 
the invention; Fig. 27 is a schematic . diagram that shows 
changes of normal deviation of distribution in the z-axis 
direction upon making a neural network by introducing an 
anisotropy merely in the z-axis direction on a 
three-dimensional lattice according to the invention; Fig. 
28 is a schematic diagram that shows changes in number of 
adjacent sites among fractal figures for a growth step during 
fabrication of a neural network by introducing an anisotropy 
merely in the z-axis direction on a three-dimensional lattice 
according to the invention; Fig. 29 is a schematic diagram 
that shows a result of simulation for growth of a coupled 
fractal network according to the invention; Fig. 30 is a 
schematic diagram that shows a result of simulation for 
growth of a coupled fractal network according to the 
invention ; Fig . 3 1 is a schematic diagram that shows a result 
of simulation for growth of a coupled fractal network 
according to the invention; Fig. 32 is a schematic diagram 
that shows a result of simulation for growth of a coupled 
fractal network according to the invention; Fig. 33 is a 
schematic diagram that shows a result of simulation for 
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growth of a coupled fractal network according to the 
invention; Fig. 34 is a schematic diagram that shows a result 
of simulation for growth of a coupled fractal network 
according to the invention; Fig. 35 is a schematic diagram 
that shows changes in distance between centers of gravity 
for a growth step in simulation of growth of a coupled fractal 
network according to the invention; Fig. 36 is a schematic 
diagram that shows changes in distance between centers of 
gravity for a growth step in simulation of growth of a coupled 
fractal network according to the invention; Fig. 37 is a 
schematic diagram that shows changes in distance between 
centers of gravity for a growth step in simulation of growth 
of a coupled fractal network according to the invention; 
Fig. 38 is a schematic diagram that shows changes in distance 
between centers of gravity for a growth step in simulation 
of growth of a coupled fractal network according to the 
invention; Fig. 39 is a schematic diagram that shows changes 
in distance between centers of gravity for a growth step 
in simulation of growth of a coupled fractal network 
according to the invention ; and Fig . 40 is a schematic diagram 
that shows changes in distance between centers of gravity 
for a growth step in simulation of growth of a coupled fractal 
network according to the invention. 

Best Mode for Carrying Out the Invention 

Explained below are embodiments of the invention 
with reference to the drawings , 



(1) Dielectric breakdown model 

First explained is a method of making a singular 
fractal figure (for example, a figure such as tree-like 
projections) ((1) A. Erzan, L. Pietronero, A. Vespignani, 
Rev. Mod. Phys. 67, 545(1995); (2) T. A. Witten, Jr. andL.M. 
Sander, Phys. Rev. Lett. 47, 1400(1984); Phys. Rev. B 27, 
5686(1983)). This method is a dielectric breakdown model 
proposed by Niemeyer et al . ( ( 3 ) L . Niemeyer , -L i Pietronero , 
H.J. Wiesmann, Phys. Rev. Lett. 52, 1033(1984)). 

As an example, a tetragonal lattice S in a 
two-dimensional space is defined, and a, scalar potential 
field <p(i, j) is defined on a lattice site (ii, i 2 , i3)eS. 
This is called an electric potential. Let this electric 
potential follow a Laplace equation: 

A^(z\i) = 0 (1) 

The figure T n defined therefrom is a set of lattice sites 
on a two-dimensional lattice. T 0 consists of (0, 0) alone, 
and T n+ i is created by sequentially adding one lattice site 
to T n by the following rule . 

Let the electric potential of each site contained 
in T n be 1 and the electric potential at a site at infinity 
be 0 . That is , 

<f)(i,j) = 0 when > co (2) 
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4>{i y j) = 1 when (ij) G T n 



(3) 



Equation (1) is solved under that boundary condition, and 
electric potential of each lattice site is determined. The 
lattice site to be added to T n to form T n+i is not included 
in T n , and it is selected from a set U n of lattice sites 
closest to T n . The number of lattice sites contained in 
U n is written as N n . 

For each site contained in U n (im, jm) (where m=l, 
2, N n ), its electric field intensity is defined as: 

E^CC) = \<f>(i»Jm) -lp" < 4) 

The probability that a site (i m , j m ) in U n is proportional 
to its electric field intensity E m (a) . That is, the 
probability is: 



(5) 



By repeating the above operation, T n is formed 
progressively. An ideal fractal will be the infinitely 
repeated extreme set : 



T oa = lim T n (6) 

In case of a=l, the above coincides with the result of 
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generation of a figure by diffusion limited aggregation ( ( 2 ) 
T . A . Wit ten, Jr. and L.M. Sander, Phys. Rev. Lett. 47, 
1400(1984) ); Phys. Rev. B 27, 5686(1983)). 
(2) Fractals that grow while interacting with each other 
Defined below are fractals that grow while 
interacting with each other (interacting fractals), i.e. 
a coupled-fractal network. As an example, fractals made 
up of N c species are taken on a tetragonal .lattice S in a 
two-dimensional space. A scalar potential field <p{±, j) 
is defined on a lattice site ( i x , i 2 , ±3 ) eS , and this is called 

a potential. Then, V x) (i, 3 K i|r (2) (i, 3) iJ/ (Nc, (i, 3) 

are also defined. They satisfy the differential equations : 

&<f>(ij) = 0 (7) 

AtfW(s,j)=. 0 ( 8 ) 
Ai/> (2) (i,i)= 0 (9) 



A^(i,j)= 0 (io) 

The figure T n defined therefrom is a set of lattice sites 
on a two-dimensional lattice, and respective lattice sites 
are classified into N c species. That is, 
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Nc 

r n =UQ ( n j) t 11 * 

j=l 

and respective species are exclusive from each other. That 
is . 

consists of a single lattice site (i\S,jiS) 
exclusively, and T n+1 is created by sequentially adding one 
lattice site to T n by the following rule. First, Equation 
(7) is solved under the boundary conditions: 

^j-)z=0 when — » oo < 13 > 

<f>(ij) = l when (i, j) G T n (14) 

and potential of each lattice site is determined. Further, 
Equations (8) through (10) are solved under the boundary 
conditions : 

tj)i k) {ij) = 0 when — > oo (15) 

tpi(i,j)=l when (iJ)eQl k) U6) 
0t*)(i,j) = -l when {ij) e 'Q™ (k ? I) (17) 
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and is determined. The lattice site to be added 

to T n to form T n+ i is not included in T n , and it is selected 

from a set o f lattice sites closest to QW . 

The number of lattice sites contained in Ujf^ 

is written as • That is, the lattice site to be added 

n 

to T n is selected from the following set: 

tfn-U^ _ - (18) 

k = l 

and the number of lattice sites contained in the set, i.e. 
the number of candidates , is : 

^ = X>< fc > (19) 

Jfc=l 

For each site (i^J 1 ^) (m-1. 2 N'l k) \ 

contained in £/(*) , intensity of its electric field 
is defined as: 

E%Ha,0) = - ir x i* , * ) (ii.".iif ) ) - n* < 20 > 

The probability that a site (imKj^) in U n is 
selected is proportional to intensity of its electric field 
£^(0,/?) . That is. the probability is: 

jfiW) (2D 

Jb=l i=i 
14 



By repeating the above-mentioned procedures, T n is formed 
progressively . 

Here is given a physical (or physiologic) 
interpretation regarding the above model . cp ( i , j ) provides 
the probability of arrival of a growth material S rc 
transported from a remote site by diffusion at a region 
permitting growth from interpretation of the dielectric 
breakdown model. On the other hand, iM±, 3) provides the 
probability of arrival of any substance X transported to 
a region permitting growth of the k-th species from portions 
where species except for the k-th species have been already 
grown. The assumption that the probability that growth of 
the k-th species occurs is proportional to the product means 
the assumption that the growth occurs only when both the 
growth material S rc and the substance X have arrived. That 
is, here is assumed the situation in which growth occurs 
by adhesion of the growth material S rc to the portion where 
the k-th species has already grown with the aid of an 
adhesive substance X. Of course, it means a reaction that 
is rate-determined by adhesion interposing the adhesive 
substance X. In an alternative interpretation, it is 
possible to consider that cf>(i, j) provides the probability 
of arrival of the growth material S rc transported by diffusion 
from a remote site to the region where it can grow, similarly 
to the above interpretation whereas \Mi, J) provided the 
potential at the position of the k-th species, and the 
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probability that growth of the k-th species occurs is 
proportional to the product of the probability of arrival 
of the growth material S rc and the intensity of the electric 
field. 

Explained below is a specific example of the method 
of fabricating a neural network based on the model shown 
in (2) above. One of results of simulation is shown in Fig. 

1, In this simulation, (a, 3) =(0.5, l.QKwas used, and 
4000-step growth from four sites was conducted. In Fig. 

2, using the same (a, p) , 3000-step growth from three sites 
was conducted. By changing parameters , various figures can 
be created. If a corresponding to the growth material 
transported from a remote site is large, then respective 
elements tend to spread. If 3 corresponding to the adhesive 
substance nearby is large, the coupling becomes dense. 
Figures upon executing 2000 -step growth from two sites by 
changing (a, 3) are shown in Figs. 3 through 7. In Fig. 

3, similarly to the above two examples, (a, 3) =(0.5, 1.0) 
was used. Examples increasing 3 from that are shown in Figs . 
4 through 6. (a, 3)=(0.5, 1.5) was used in Fig. 4, (a, 
3)=(0.5, 2.0) was used in Fig. 5, and (a, (3) = (0.5, 2.5) was 
used in Fig. 6, respectively. It is appreciated that the 
coupling becomes dense as 3 increases. On the other hand, 
if the relation in value between a and 3 is reversed, growth 
tends to proceed spreading remoter as shown in Fig. 7 [(a, 
3) = (1.0, 0.5)]. In this manner , by changing the parameters 
a and 3, various neural network structures can be produced 
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while controlling properties of figures. 
(3) Extended interacting fractals 

Extended interacting fractals, i.e. a 
coupled- fractal network, are defined below. As an example , 
fractals made up of N c species are taken on a tetragonal 
lattice S in a two-dimensional space. A scalar potential 
field (J>(i, j) is defined on a lattice site (i, j)€S, and 
this is called a potential. Then, i|r (1) (i, j ) , V 2) (i, 3)* 
.... 4f (Nc) (i, 3) are also defined. They satisfy the 
differential equations: 

A<f>(iJ) = 0 (23) 
A^< 1 >(i,j)= 0 (24) 
A^( 2 >(i,j)= 0 (25) 



A%l>( N ^(iJ)= 0 (26) 

The figure T n defined therefrom is a set of lattice sites 
on a two-dimensional lattice, and respective lattice sites 
are classified into N c species. That is, 

T n =\jQP < 27 > 
i=i 
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and respective species are exclusive from each other. That 
is. 



■Qti } r\QP = *-*j¥>k (28) 

consists of a single lattice site (t^ih^) 
exclusively, and T n+i is created by sequentially adding one 
lattice site to T n by the following rule. First, Equation 
(23) is solved under the boundary conditions: 

<f>(ij) = 0 when (i, j) — > oo (2g) 

<f>(ij) = l when (hj)€T n (30) 

and potential of each lattice site is determined- Further, 
Equations (24) through (26) are solved under the boundary 
conditions: 

^ k) (iJ) =^oo when (i, j) > oo (31) 

= l when (hj) e Qn ] (32) 

^ k) (iJ) = -1 when (ij) e Q« (k ±1) (33) 

and is determined. The lattice site to be added 

to T n to form T n+i is not included in T n , and it is selected 
from a set of lattice sites closest to Qffl • 
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The number of lattice sites contained in 
is written as • That is, the lattice site to be added 

to T n is selected from the following set: 

U n =\J U i k) (34) 
fc=l 

and the number of lattice sites contained in the set, i.e. 
the number of candidates , is : 

Nn=T r ^ ) OS) 
k=l 

For each site contained ir 

intensity of its electric field is defined as: 

= MWJP) - If >< W'HWJP) ~ IP (36) 

The probability that (i^iim*) in u " is 
selected is proportional to intensity of its electric field 
EW(oc,/3) . That is. the probability is: 



p£W) = X^W) (37) 



A 

A = EE^'M (38) 

By repeating the above-mentioned procedures, T n is formed 
progressively . 

In the above model, the model explained in (2) is 
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expanded in the portion of the boundary condition 
ipW(i,j) = ^oo o f ^ k Khj) ' That is ' by limitation of 
11/00=0, the above model results in the model of (2). 

Explained below is a specific example of the method 
for fabricating a neural network based on the model in ( 3 ) 
above. Results of simulations are shown in Figs, 8 through 
14. In these simulations, (a, P)=(0.5, 1.0) was used, and 
3000 -step growth from three sites was conducted. Fig. 8 
shows a case of growth under the boundary condition of Hf~=0, 
and it is the same as (2) above. i|f.= -l was used in Fig. 
9, \M=-0.6 was used in Fig. 10, and i|/oo=-0 . 2 was used in Fig. 
11. 

In case that ilr« is negative, the potential 
difference between fractal figures belonging to species 
other than itself and a site at infinity becomes smaller, 
and the interaction among fractal figures becomes weaker . 
On the other hand. Figs. 12 through 14 show those in which 
Uf« is a positive value. ^«=0 was used in Fig. 12, \M=0.6 
was used in Fig. 13, and \|/«,=0.2 was used in Fig. 14. 

In case of i|f«>0 . in contrast, interaction among 
fractal figures becomes stronger, and a shape massed more 
densely is obtained. 

(4) Extended interacting fractals introducing anisotropy 
Extended interacting fractals introducing 
anisotropy, i.e. a coupled-fractal network, are defined 
below. As an example, fractals made up of N c species 
are taken on a tetragonal lattice S in a two-dimensional 
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space. Regarding S as an anisotropic field, an anisotropic 
parameter tensor M=(m^, v ) , L=(l u . v ) is introduced here . A 
scalar potential field $ M (i, j ) is defined on a lattice site 
(i, j)eS, and this is called a potential. Then, 

^L\iJ)^L\hj)^-'^L C) ^j) ^ alS ° d6fined - 

They satisfy the differential equations: 

M***™^** ' ■ - - (6i) 

■ ; ^,v dr » or u 

- 91 ° tPfrj) = 0 (63) 



u f v=x,y 3 r t* 



M)l /=x,y 'm," c/7 >' 

Here , represents a difference on the lattice site, and 



d_ 

for example. 



d<j> M {i,j) _ 4>Af(i + I J) - <t>M(iJ) 

dr x ~ (* + !)-« (65) 



The figure T n defined therefrom is a set of lattice sites 
on a two-dimensional lattice, and respective lattice sites 
25 are classified into N c species. That is, 

T n =\jQP (66) 

p=l 
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and respective species are exclusive from each other, 
is , 



That 



HOP = *P*<1. (67) 



Qq^ consists of a single lattice site 

(Ap) Ap)\ exclusively, and T n+i is created by sequentially 
\*ini > Jim ) 

adding one lattice site to T n by the f ollowing xule . First , 
Equation (61) is solved under the boundary conditions: 

<£ M (z,j) = 0 when (ij) > oo (68) 

<f>M(iJ) = l when (iJ)eT n (69) 

and potential of each lattice site is determined . Further , 
Equations (62) through (64) are solved under the boundary 
conditions : 

$l\*J) = when (ij) >oo ( 70 } 

V>i p) (z,j) = l when (i,j)€g?>. (71) 

*l\*J) = -i when (ij) e (P ± q) < 72 ) 

and j) is determined. The lattice site to be added 

to T n to form T n+ i is not included in T n , and it is selected 
from a set of lattice sites closest to Q(p) . 
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The number of lattice sites contained in f/W 
is written as . That is, the lattice site to be added 

to T n is selected from the following set: 

Un=\JUjp (73) 
P=l 

and the number of lattice sites contained in the set. i.e. 
the number of candidates . is : 

N n = Yl N n ) < 74) 
P=l 

For each site (m=l, 2. . . . , jV^) ) 

contained in Ujjfi , intensity of its electric field is 
defined as : 

i#>(a,/?,M,Z) - - ir x |^>-(iW,iW) - If 

(75) 

The probability that a site in U„ is 

selected is proportional to intensity of its electric field 
E$(ct,P,M,L) . That is, the probability is: 



P£>(<*,/?,M,L) = li#>(a,/?,M,£) 



(76) 



N ' N » P) (77) 
p=l m=l 

By repeating the above-mentioned procedures, T n is formed 
progressively . 
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The above model has been extended from (3) above 
as a result of introduction of anisotropic parameters in 
Equations (61) through (64). By limitation of 
m /£ , v =m5 )£i v , l fttV ~15 Ht v (m, l=const . ) , the above model results 
in the model of (3). 

Explained below is a specific example of the method 
for fabricating a neural network based on the model in (4) 
above , Results of simulations are shown below. _ First shown 
is how the coupled-fractal changed. For easier 
understanding, anisotropy is introduced in only one 
direction on a two-dimensional lattice (in this case, in 
the y-axis direction on an x-y plane. That is. 

While fixing the parameters as (a, p)«(0.5', 1.0) and \J/„=0, 
4000-step growth was conducted from two sites on the 
y-axis. Figs. 15 through 20 are coupled-fractal figures 
upon m=l=0 . 5 , m=l=0 . 75 , m=l=1.0, m=l=1.25, m=l=l . 5 and 
m=l=2, respectively. It is appreciated from these figures 
that as m and 1 increase, the coupled-fractal figure tends 
to extend in the y-axis direction (growth in the y-axis 
direction is promoted) and that as m and 1 decrease, the 
coupled-fractal figure tends to extend in the x-axis 
direction (growth in the y-axis direction is restricted). 
It is further known from these structural changes that the 
number of adjacent sites among different fractal figures 
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changes. More specifically, in the figure shown in Figs. 
15 through 17, there are no adjacent sites between every 
two fractal figures. However, adjacent sites gradually 
increase from six in Fig. 18, through 57 in Fig. 19 to 95 
in Fig. 20, for example. In this manner, by changing 
anisotropic parameters , the coupling mode among a plurality 
of fractal figures can be changed. Fig. 21 shows a figure 
as a result of growth at different values of m=0 . 5 and 1=2 . 0 . 
It is appreciated from it that influences of the anisotropic 
parameter 1 introduced into are larger than influences 
of the anisotropic parameter m introduced into <$> . As readily 
understood from the fact that the coupling strength changes , 
the method by this model is not mere spatial expansion and 
contraction, but it is an improvement by anisotropy of growth 
algorithm introduced absolutely first. 

(5) Extended interacting fractals introducing anisotropy 
in a three-dimensional space 

Extended interacting fractals, i.e. a 
coupled-fractal network, developed from the model in (4) 
above and introducing anisotropy in a three-dimensional 
space is introduced. Fractals made up of N c species are 
taken on a tetragonal lattice S in a two-dimensional space. 
Regarding S as an anisotropic field, an anisotropic 
parameter tensor: 
M={m u , v ) , L=(l ; ,, y ) 

is introduced here. A scalar potential field (J>M(i# 3* ^) 
is defined on a lattice site (i, j, k)es, and this is called 
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a potential. Then. *£ ) «.j. *) , ^\iJ, *) , • • -, ^(i.i.*) 

are also defined. They satisfy the differential equations : 



ir — —SrMhj^) - 0 



(79) 



ix,i/=x,t/ ; z ar p hh^ ° r v 

d 1 a , (2) 



fi,v=x,y,z d r f* ® r *> 



*I>L (v**). = 0 (81) 



d 

Here, Q r represents a difference on the lattice site, 
and for example, 

d<f>M(iJ, k) _ 4> M {i + *-) ~ Mi J, {83) 
dr x (*.+ !)-* 

The figure T n defined therefrom is a set of lattice sites 
on a two-dimensional lattice, and respective lattice sites 
are classified into N c species. That is, 

N c 

T n = (JQ^ (84) 
P=l 

and respective species are exclusive from each other. That 
is , 



gl p) n^ )==0 ' ifp** (85) 
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Vo consists of a single lattice site 
(4vlyJMj kin!) exclusivel Y * and T n+i is created by sequentially 
adding one lattice site to T n by the following rule. First, 
Equation (79) is solved under the boundary conditions: 

<£ M (i,j,fc) = 0 when (ij,k) — ► oo (8 6) 

<f>M{iJ,k) = l when (i>j>k)€T n (8 7) 

and potential of each lattice site is determined. Further, . 
Equations (80) through (82) are solved under the boundary 
conditions: 

^£\^h k ) = ^oo when (ij,k) >oo (88) 

tfi p) (*Vi,*) = 1 when (i.j.Jfe) <E Q« (89) 
0l p) («,i,*) = -l when *)€<?£?(?#«) (90) 

! 

and 'f/^\ijjjk) is determined. The lattice site to be added 
to T n to form T n+i is not included in T n , and it is selected 
from a set Jj(p) of lattice sites closest to Q(p) . 

The number of lattice sites contained in 7/00 

Tl 

is written as ]V(p) . That is, the lattice site to be added 

n 

to T n is selected from the following set: 
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u n = \jup (91) 

and the number of lattice sites contained in the set, i.e. 
the number of candidates , is : 

^n=E^ P) (92) 

For each site (i£\j£\kM) (m=l, 2, _ N^) 

contained in , intensity of its electric field is 

defined as: 

EM(«,frM,L) = |M4\#>,*M) - 1| 0 x - If (93) 

The probability that a site (ifi?,^, in U„ is 

selected is proportional to intensity of its electric field 
£j£>(a,£,M,Z,) - That is. the probability is: 

P%Ha,0 y M,L) = ±EM(<*,P,M,L) (94) 
A = EE f3,M,L) < 95 » 

P=l 771 = 1 

By repeating the above-mentioned procedures, T n is formed 
progressively. 

Explained below is a specific example of the method 
for fabricating a neural network based on the model in ( 5 ) 
above. Results of simulations on a three-dimensional 
lattice are shown below. Introducing an anisotropy merely 
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in the z-axis direction on a three-dimensional lattice, that 
is, determining as: 



and fixing parameters as (a, 3)= (0.4, 0.8) andi|/oo=0, 5000-step 
growth was conducted from three sites in an x-y plane . Figs . 
22 through 26 show coupled-fractal figures upon m=l=0.02, 
m=l=0 . 1 , m=l=l . 0 , m=l=2.0 and m=l=10 , respectively. Fig. 
27 shows a normal deviation s in the distribution along the 
z-axis direction. In Fig. 27, the abscissa shows values 
of m=l=x. It is appreciated from these figures that 
distribution in the direction introducing the anisotropy 
can be controlled on a three-dimensional lattice. That is , 
by decreasing the value of m=l=x, growth in the z-axis 
direction can be suppressed, and to the contrary, by 
increasing the value of m=l=x, growth in the z-axis direction 
can be promoted . Here is reviewed that the number of ad j acent 
sites among different fractal figures can be controlled 
approximately. Fig. 28 is a diagram made by plotting the 
number N of adjacent sites under the same situation as Figs . 
22 through 26, namely, under m=l=0 . 02 , m=l=0.1, m=l=1.0, 
m=l=2.0 and m=l=10, taking growth steps on the abscissa. 
It is appreciated from Fig . 28 that as values of the parameters 
m and 1 decrease, i.e. , as the growth sites extend in the 
x-y plane defining the start sites, the number of coupling 
progressively increases . 



M = 



( 1 ° 
0 1 

^ 0 0 




(96) 
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The simulation program for fabricating the neural 
network can be supplied in form of a recording medium that 
can be read through a computer. 

(6) Statistical interaction in a coupled-fractal network 
It is natural that fluctuation occurs in a structure 
made through a random process or in the processing using 
the structure due to the random process . The same applies 
to information processing in our brains. However, 
according to the matters having been clarified along with 
recent developments of the statistical physics, there is 
a phenomenon or a universal nature that appears in such a 
random system only when its statistic distribution is 
averaged . 

In a coupled fractal network, if a row of random 
numbers is changed in its growth process, then the shape 
itself of the network to be made will change, but universal 
properties independent from the difference are believed to 
exist. One of them is the statistical interaction defined 
by the statistic distribution. The nature there of is 
clarified here. 

Growth of a coupled-fractal network is 
characterized by two parameters (a, 3) . However, in actual 
growth tests, various networks are formed due to random 
variables. A review is made for universal properties not 
affected by the random variables. 

In the simulation made below, growth of M c =2 was 
conducted on a two-dimensional tetragonal lattice of 
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201x201. The first start site that is the element of 

Qq^ is (185, 100), and the second start site that is the 
element of is (117, 100). Regarding growth of n=800 

steps, growth simulation was conducted for a different row 
of random numbers of M=300. 

Examples of simulation results are shown in Figs. 
29 through 34. (a, p)=(0.4, 0.6) was used in Fig. 29, (a, 
3)=(0.8, 0.6) was used in Fig. 30, (a, (3) = (1.2, 0.6) was 
used in Fig. 31, (a, 0)=(O.4, 1) was used in Fig. 32, (a, 
(3) = (0.8, 1) was used in Fig. 33, and (a, p)=(1.2, 1) was 
used in Fig. 34, respectively. In each of these figures, 
four examples of coupled- fractal networks are shown. They 
were grown by using the same (a, (3) and different rows of 
random numbers. Every four examples are different 
coupled-fractal network structures, but have common 
properties, that is, properties regulated by (a, p). One 
of them is the fractal property of each element fractal, 
and it does not rely on interaction. What is discussed here 
is interaction among element fractals. 

Interaction among element fractals can be defined 
solely by using a statistic average concerning a plurality 
of growth executed by using different rows of random numbers 
because there is no way of discriminating from a single growth 
simulation result whether the relationship among element 
fractals is regulated by interaction or individual element 
fractals have been formed independently Therefore, 
interaction among element fractals occurring in the 



statistic average of a number of growth simulations is called 
statistical interaction. 

Growth of M=300 was conducted, and let each growth 
be distinguished from others by writing it with the suffix 
p as T n (p), p=l, 2, , . -M. It can be written as: 

Tn(p) = U QPiP) (lOD 

by way of the element Q^(p) of each coupled- fractal 
network. The number of elements of i.e., the 

number of sites where k-th species in the p-th growth 
simulation are grown in n steps is written here as M n ,k. P - 
Let the barycentric coordinates of Q^(p) be introduced 
into : 



W n* } (p) = TT S r (102) 



The sample average of barycentric coordinates is : 

! M 

w<*> = -£I><*>(i>) ' u°3> 

iV1 P =l 

and the distance between the gravity expectation values: 

is a convenient quantity. Then, for the purpose of analyzing 
the statistic interaction, a correlation function is 
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introduced. 



It is: 



n 



M 



M 



(105) 



p=l 



where 



1 



5D ( r " w 




(106) 



What is analyzed here is an example of N c =2 . In particular, 
here are computed the expectation value of the barycentric 
distance of first and second species, namely: 



and the dimensionless quantity as the correlation intensity 
of the first species and the second species, namely: 



Barycentric distance in case of (3 = 0.6 is shown in Fig. 35, 
that in case of (3 = 0.8 is shown in Fig. 36, and that in case 
of 3 = 1 is shown in fig. 37, respectively. It is appreciated 
that the barycentric distances increase with the increase 
of a , respectively . That is , since the barycentric distance 
always increases, it means that the coupled-fractal network 
grows outward. In this case, it is permitted to consider 
that a typical repulsive interaction is working. On the 



D(n) = i#; 2 > 



(107) 




(108) 



33 



other hand, in case that D(n) decreases with n, a condensed 
network is formed. In this case , it is permitted to consider 
that a strong attractive interaction is working. 

Here is reviewed the interaction intensity X (n) . 
Fig. 38 shows that in case of 3 = 0. 6 , Fig. 39 shows that in 
case of 3 = 0.8, and Fig. 40 shows that in case of 3 = 1- In 
the region of n<100, a large difference is not observed. 
This is the region where growth progresses independently, 
and statistical interaction among species does not affect 
the configuration. In the region of 100<n<200, X (n) 
maintains an approximately constant value around 0.2. 
Thereafter, in the region of n>200 , statistical interaction 
among species depending upon (a, 3) appears. It is 
appreciated that x decreases with an increase of a. This 
can be understood more easily by giving a physical 
interpretation of X (n). For example, let the specie 1 
deviate to the left. Then, if the species 2 is also pulled 
to the left, it results in X (n)>0. That is , if X (n) is 
a positive, large value, the species 2 is pulled by the species 

1 and follows it. On the other hand, if x(n)<0, the species 

2 results in growing in the opposite direction, and it can 
be said to be repulsive . Therefore , the interpretation that 
the attractive force gets weaker with an increase of a can 
be concluded to be quantitatively accurate in a statistical 
sense. It is appreciated that the statistical interaction 
is affected also by a difference of (3. In Fig . 40 , x changes 
largely in response to the change from a-0.2 to a=1.2, but 
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this change decreases with the decrease of (3 . In the growth 
model of the coupled-fractal network, 3 controlled the 
adsorption probability of the adhesive substance. It is 
presumed that a decrease of (3 weakens the dependency of the 
adhesive substance upon its location, and thereby decreases 
the influences from positions where the other species are 
formed. 

Summarizing the foregoing discussion, in a 
coupled- fractal network, statistical interaction among 
element fractals is defined by an ensemble average of a 
plurality of growth tests, and the universal nature of the 
network is controlled by the statistical interaction. 

Although the invention has been explained by way 
of specific examples , the invention is not limited to those 
embodiments , but can be changed or modified in various modes 
based not departing from the technical concept of the 
invention. 

As described above, according to the invention, 
since fractal structures are grown from a plurality of start 
sites, respectively, while having them interact with each 
other, to form fractal structures coupled to each other, 
a complicated network such as neural network can be made 
easily in a well-controlled manner. Especially when the 
space for growing fractal structures has an anisotropy, the 
coupling mode among different fractal structures can be 
controlled, and therefore, complicated networks such as 
neural networks, having a more variety of structures, can 
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be made easily in a well-controlled manner. 
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